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Abstract 

We discuss new physical phenomena expected in particle production in hadron-hadron collisions 
at high energy, as a consequence of Pomeron loop effects in the evolution equations for the Color 
Glass Condensate. We focus on gluon production in asymmetric, 'dilute-dense', collisions : a 
dilute projectile scatters off a dense hadronic target, whose gluon distribution is highly evolved. 
This situation is representative for particle production in proton-proton collisions at forward 
rapidities (say, at LHC) and admits a dipole factorization similar to that of deep inelastic 
scattering (DIS). We show that at sufficiently large forward rapidities, where the Pomeron 
loop effects become important in the evolution of the target wavefunction, gluon production 
is dominated by 'black spots' (saturated gluon configurations) up to very large values of the 
transverse momentum, well above the average saturation momentum in the target. In this regime, 
the produced gluon spectrum exhibits diffusive scaling, so like DIS at sufficiently high energy. 
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1 Introduction 



The recent data on particle production at semi-hard transverse momenta {p± = 1 -j- 5 
GeV) in deuteron-gold (d+Au) colhsions at RHIC [1-3] show strong evidence in favour 
of high gluon density effects in the nuclear target and can be naturally accommodated 
[4-14] within the same theoretical framework that has been previously used to describe 
the HERA data at small Bjorken-x [15] : the effective theory for non-linear evolution 
in QCD at high energy known as the "Color Glass Condensate" (CGC) [16, 17] (see the 
review papers [18] for details and more references). 

In particular, the results at RHIC show a striking difference between the particle yield 
at central and, respectively, forward rapidities — namely, a very rapid evolution of the 
nuclear modification factor i?dAu from a 'Cronin peak' (i?dAu > 1) at central rapidity 
(77 = 0) to 'high-p^ suppression' (i?dAu < 1) at forward rapidity {r] > 1); see, e.g., the 
discussion in [19] — which is naturally attributed to the small-x evolution of the gluon 
distribution in the nucleus in the presence of non-linear effects [4, 6, 8, 10, 13]. (One should 
recall at this point that increasing 1] is tantamount to probing gluons with smaller values 
of X in the nuclear wavefunction; see the discussion of kinematics in Sect. 2.) Moreover, 
for very forward rapidities (77 > 3), the spectrum of the produced hadrons appears to be 
consistent with geometric scaling [20], which is a hallmark of the BFKL evolution in the 
vicinity of the saturation line [21-23]. For smaller rapidities < r] < 3, the violations of 
this scaling become visible too [14], and they indeed follow the pattern predicted by the 
CGC effective theory [21-23] — a pattern which has been also verified in the HERA data 
for the F2 structure function at a; < 0.01 [15]. 

These results at RHIC demonstrate the potential of the 'dilute-dense' (here, deuteron- 
gold) hadron-hadron collisions to accurately probe the universal phase of QCD at high 
energy — the Color Glass Condensate. We expect this potential to be amplified in the 
incoming experiments at LHC, because of the considerably higher energies that will be 
available there. One should emphasize here that such asymmetric collisions, in which one 
of the partaking hadrons (the 'projectile') is dilute and the other one (the 'target') is 
dense, are exceptional in several respects: 

• On the experimental side, the 'dilute-dense' collisions are less affected by the final-state 
interactions than the 'dense-dense' collisions (like the nucleus-nucleus ones), hence they 
represent a relatively clean measurement of the initial conditions (nearly on the same 
footing as the deep inelastic lepton-hadron scattering). 

• On the side of the theory, such collisions are better under control, in the sense that 
factorization schemes are available, which allow one to compute the cross-section for 
particle production as a convohition of (standard) parton distributions for the dilute 
projectile times 'generalized parton distributions' (including high-density effects) for the 
dense target times partonic cross-sections [9,24-29]. This factorization is very similar to 
the 'dipole factorization' used for DIS at high energy (see, e.g., Ref. [18] and references 
therein), and the 'generalized parton distributions' alluded to above are in fact dipole- 
hadron scattering amplitudes [30] , which describe the multiple scattering between a 'color 
dipole' and the dense hadronic target. 
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Hence, the non-trivial scattering problem that one has to solve is that of the scattering 
between an elementary color dipole and the high-density hadronic target. To that aim, 
one can rely on the CGC formalism [18], where the dipole amplitudes obey an infinite 
hierarchy of non-linear evolution equations — the Balitsky-JIMWLK hierarchy [17, 31- 
33] . In the limit where the number of colors Nc is large, this hierarchy boils down to a single 
equation, the Balitsky-Kovchegov (BK) equation [31,34]. Most of the previous analyses 
of the d+Au collisions within the framework of CGC were based on the solution to BK 
equation, or on various approximations to it [4, 6, 8, 10]; as mentioned before, they led to a 
rather satisfactory picture, which provides a qualitative explanation^ of the experimental 
situation at RHIC and is robust against the various approximations. 

But whereas the current CGC formahsm, centered around the JIMWLK equation and its 
main consequences like 'geometric scaling', appears to be appropriate for understanding 
the experimental situation at RHIC and HERA, it is not clear whether this success will 
extend to the future data at LHC. Indeed, recent theoretical developments [35-45] show 
that some physical ingredients which have been ignored by the JIMWLK evolution play 
an increasingly important role with increasing energy, and should lead to a dramatically 
different physical picture at sufficiently high energy. These ingredients are the gluon- 
number fluctuations [35, 37, 46, 47] associated with gluon splitting (bremsstrahlung) in 
the dilute regime, which act as a seed for K-body correlations with k > 2. Although 
formally suppressed by higher powers of the coupling constant as, such correlations are 
rapidly amplified by the BFKL evolution (the faster the larger is k) and eventually play 
an essential role in the evolution towards gluon saturation at high energy. The Balitsky- 
JIMWLK equations correctly encode the role of such many-body correlations in providing 
saturation, but they fail to include the physical source for such correlations, namely, the 
gluon-number fluctuations at low density [38]. Thus, while these equations are properly 
describing the evolution of a dense system towards saturation, they cannot describe the 
formation of such a dense system via evolution from a dilute system at low energy. 

At this point, one may be tempted to conclude that the Balitsky-JIMWLK equations 
should at least apply to the evolution of a large nucleus, which starts with a relatively 
high gluon density already at low energy. However, this is not true either ! Even for a 
large nucleus, the gluon distribution involves a dilute tail at relatively high transverse 
momenta (well above the nuclear saturation momentum), and the evolution of that tail is 
still dominated by fluctuations. It is nevertheless true that, in order to probe the effects of 
fluctuations for such a large nucleus, one needs to go up to much higher energies, namely, 
high enough for the saturation effects to have propagated at relatively large transverse 
momenta, which were originally in the dilute tail at low energy. In that sense, the evolution 
of a originally dilute system, like a proton or a color dipole, is a better laboratory to study 
the effects of ffuctuations. 

The equations describing gluon evolution in the presence of both saturation and fluc- 
tuations are presently known only for large and are generally referred to as the 

^ Quantitative fits have been also obtained [11, 13, 14], but only at the expense of introducing 
some free parameters, to account for the lack of accuracy of the current formalism and for the 
uncertainties associated with the non-perturbative initial conditions. It is interesting to note 
that most of these parameters were in fact fixed from fits to the HERA data for F2 . 
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'Pomeron loop equations' [38-40,42]. The 'Pomeron' here is the BFKL Pomeron [48]; 
namely, this is the amplitude for the scattering between an external dipole and the tar- 
get in the linear, leading-logarithmic, approximation at high energy. In an appropriate 
gauge, this amplitude has a diagrammatic interpretation in terms of ladder diagrams 
whose rungs are strongly ordered in the 'rapidity' variable Y = ln(l/a:;). At large Nc, 
the gluon bremsstrahlung (responsible for gluon number fluctuations) can be described as 
1 — s> 2 Pomeron splitting, while the saturation effects are tantamount to 2 ^ I Pomeron 
merging. Hence, the equations describing the complete evolution involve both (Pomeron) 
splitting and merging, and thus generate Pomeron loops through iterations. 

The Pomeron loop equations form an infinite hierarchy, with a complicated non-local 
structure, and in spite of intense theoretical efforts, their general solutions are not yet 
known. Still, the asymptotic behaviour of the solutions at high energy is known because 
of the universality of the stochastic process described by these equations [37, 49-51] : this is 
in the same universality class as the 'reaction-diffusion' process A ^ 2A, which is relevant 
to a large variety of problems (in statistical physics, chemistry, biology ...) and has been 
intensely scrutinized over the last decades (cite Ref. [52] for a review). But even in that 
context, it was only recently realized that this process is extremely sensitive to particle- 
number fluctuations in the dilute regime [53]. In the language of QCD, this sensitivity 
explains why the predictions of the Pomeron loops equations turn out to be very different 
from those of the Balitsky-JIMWLK equations, which ignore fluctuations [36-38]. 

By using known results from statistical physics [52, 53], it has been possible to deduce the 
functional form of the dipole scattering amplitudes in QCD at high energy and large N^, in 
terms of juste a couple of uncontrolled parameters [38]. The most striking feature of these 
results is the fact that, even in the weak scattering regime, that is, for a relatively small 
projectile dipole, the scattering amplitude is dominated by black spots, i.e., rare gluon 
conflgurations which are at saturation on the resolution scale ~ of the projectile 
(r is the dipole size) and thus look 'black' to the latter — the dipole is completely absorbed 
when hitting any such a spot. The reason why the average amplitude is nevertheless small 
is because such 'black spots' are relatively rare, so the target disk looks transparent (or 
'white') on that resolution scale at most impact parameters (see Sect. 3 for details). 

Very recently, the consequences of this new physical picture have been explicitly worked 
out for the case of DIS [45] , with the conclusion that new phenomena are to be expected 
at very high energy — chiefly among them, the replacement of geometric scaling by a new, 
diffusive, scaling which should extend up to very large (cf. Sect. 3). But whereas in 
the case of DIS, it is not clear when, and whether, new experiments will become available 
with sufficiently high energy to probe this new physics (as aforementioned, the current 
experimental situation at HERA is consistent with geometric scaling, thus suggesting an 
intermediate-energy regime), the imminent advent of the LHC opens the possibility to 
observe the Pomeron loops via high-energy hadron-hadron collisions. 

Both pp and pA (and also AA) collisions will be performed at LHC, but for the rea- 
sons explained earlier, the best process to study the fluctuation-dominated regime at 
high energy is particle production at forward rapidities in pp collisions. Then, one of the 
participating protons (the 'projectile') is dilute and acts as a probe of the small-x part 
of the wavefunction of the other proton (the 'target'), which is highly evolved. Ideally, 



4 



one needs very large forward rapidities for the produced particles, in order to generate a 
large asymmetry between the evolutions of the target and the projectile, and thus probe 
Pomeron loop effects in the target proton while keeping the projectile proton dilute. 

In view of the theoretical uncertainties alluded to above, it is not possible for us today to 
make quantitative predictions for LHC, and not even to reliably predict whether the new 
physics should manifest itself within the kinematical range covered by the LHC, or not. 
What we can do, however, is to predict new, qualitative, phenomena which would unam- 
biguously signal the new physics in the data for particle production, and even compute 
these phenomena in terms of few free parameters. By comparing these predictions to the 
data to become later available at LHC, one should be able to identify this new physics (if 
it is present indeed !) and extract the values of the relevant parameters. 

Our main purpose here is to demonstrate the consequences of the Pomeron loop physics 
on the particle production in hadron-hadron collisions at high energy. To that aim, we 
shall focus on a problem which is conceptually transparent and avoids unnecessary com- 
plications ^ : that of the forward gluon production in onium — hadron scattering at high 
energy and large Nc- 

• The 'onium ' is a dilute hadronic system produced via the BFKL evohition of a small 
color dipole. Its wavefunction is explicitly known within perturbative QCD at large Nc 
[46], and this is why we have chosen it as a projectile. But the generalization of our 
subsequent results to a more realistic projectile, like a (dilute) proton, is straightforward: 
In the most interesting kinematical region, the projectile is merely represented by the 
standard (integrated) gluon distribution function, as determined by the leading-twist, 
DGLAP, approximation. 

• The hadronic target will be assumed to be probed at very small values of x (corre- 
sponding to very forward rapidities for the produced gluon), in the fluctuation-dominated 
regime which is universal: this is the ultimate form of hadronic matter which is produced 
at sufficiently high energy starting with arbitrary initial conditions at low energy. In this 
regime, the dipole-target scattering amplitude can be taken in the simplified functional 
form predicted by the correspondence with statistical physics, which represents an ap- 
proximate solution to the Pomeron loop equations. This functional form describes not 
only the diffusive scaling regime at asymptotically high energies, but also the geometric 
scaling regime at intermediate energies, and the transition between the two. Hence, by 
increasing the rapidity of the produced gluon, we shall be able to study the increasing 
influence of fluctuations on the cross-section for particle production and the interplay 
between geometric and diffusive scaling. 

This paper is organized as follows: The first two sections contain no new results, but 
rather collect some previous results — concerning the dipole factorization of forward gluon 
production (in Sect. 2) and, respectively, the calculation of the dipole amplitudes in the 
presence of Pomeron loop (in Sect. 3) — which are necessary for the subsequent analysis. 
The first issue (the factorization of the gluon production) is rather well estabhshed in 

^ Like the use of realistic parton distributions or fragmentation functions; note, however, that 
such 'complications' may become essential and unavoidable in view of realistic applications to 
the phenomenology. 
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the literature and will be only succinctly described here. But the second one (the high- 
energy evolution with Pomeron loops) has met with important developments over the 
last two years, whose results are perhaps less well known (since spread over several, more 
technical, papers [36-38,45]), yet they are essential for our present analysis. To cope with 
that, in Sect. 3 we shall try to give a self-contained and pedagogical discussion of these 
recent results. Then, in Sects. 4 and 5, we shall apply them to the analysis of forward 
gluon production at high energy. We shall separately discuss the relevant 'generalized 
gluon distribution' (a special Fourier transform of the dipole amplitude) (in Sect. 4) and 
the cross-section for gluon production (in Sect. 5). We shall point out some interesting 
differences with respect to the corresponding analysis of DIS in Ref. [45]. But in spite 
of these differences, the final conclusions are qualitatively similar: the cross-section for 
forward gluon production is dominated by 'black spots' up to very large values of the 
gluon transverse momentum and the ensuing spectrum shows diffusive scaling. We shall 
conclude, in Sect. 6, with a short summary and a list of open problems. 



2 Gluon production in onium— hadron scattering: Factorization 

In this section, we shall more precisely describe the calculation that we intend to perform 
and recall the factorization scheme which lies at the basis of this calculation. 

As explained in the Introduction, we are interested in gluon production in the asymmetric 
collision between a dilute projectile, whose wavefunction has been evolved in rapidity up 
to yi, and a dense target, with rapidity y2 = Y — Hi ^ Hi- {Y is the total rapidity gap 
between the projectile and the target.) Here, by "dilute" and "dense" we understand that, 
on the resolution scale set by the transverse momentum k±_ of the produced particle, the 
gluon distribution in the projectile can be described in the leading-twist approximation 
(i.e., its evolution up to yi lies within the realm of the linear, DGLAP or BFKL, evolu- 
tion equations of QCD), whereas the gluon distribution in the target is rather the site 
of important non-linear phenomena, and is strongly influenced by both saturation and 
gluon-number fluctuations — hence, by Pomeron loops. 

For definiteness, and also for more similarity with the physically interesting case of pp 
coUisions at LHC, we shall consider that the coUision is viewed in the center-of-mass 
frame, and that the two partaking hadrons are identical in that frame — or, more precisely, 
the /oo A; identical when probed via particle production at centra/ pseudo-rapidity ^ , ?7 = 0. 
However, differences between the target and the projectile will appear when one considers 
particle production at forward (pseudo)rapidity > 0, since in that case one probes rather 
different gluon evolutions in the respective wavef unctions. Speciflcally, a gluon emerging 
with rapidity r] and transverse momentum k±_ has been produced via the fusion of a 
projectile gluon having longitudinal momentum fraction Xi = {k^/^/s)c^ with a target 
gluon for which X2 = {k±/ ^)e~^ . We deduce that yi = ln(l/a;i) = Y/2 — r] and similarly 

^ The pseudo-rapidity is defined as ry = — lntan(0/2), with 9 the angle between the direction of 
the produced particle and the collision axis. Our conventions are such that 'forward rapidities', 
i.e., positive values for rj, correspond to particles emitted in the fragmentation region of the 
projectile. For a produced gluon, the pseudo-rapidity coincides with the standard rapidity. 
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Fig. 1. A schematic illustration of a Feynman graph contributing to the cross-section for gluon 
production in onium-hadron scattering in the presence ofPomeron loops. The upper part of the 
graph represents the 'onium', i.e., the modulus squared of the wavefunction of the BFKL-evolved 
dipole. Each double-line denotes a soft gluon inside the wavefunction (the gluon is represented 
as a pointlike quark-antiquark pair in a color octet state, as appropriate at large Nc). The 
produced gluon (in red) is radiated by one of the dipoles, at transverse coordinate x in the 
direct amplitude and, respectively, at y, in the complex conjugate amplitude. This effectively 
builds up a gluon-gluon dipole of size r = x — y, which multiply scatters off the strong color 
fields in the target (represented in the lower part of the figure). The evolution of the target 
wavefunction included Pomeron merging and splitting, hence Pomeron loops. 

y2 = ln(l/a;2) = Y/2 + rj, with Y = ln{s/k'j_). Thus, when both Y/2 and 1] are large 
enough, but not very different from each other, one can reach the interesting situation 
where the projectile looks dilute, while the target is dense. More precisely, we shall work 
under the assumptions that 



{(is = c^sNc/n) where the last condition [35,46] ensures that saturation effects remain 
negligible in the (relevant part of the) projectile wavefunction. In fact, the most interesting 
regime for us here is the high-energy regime at dsY > ln(l/a^) (and similarly for asV)^ 
in which the target wavefunction is sufficiently evolved to be sensitive to Pomeron loops. 



F, r/ > — 




2 ' 



(2.1) 
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Although there is no real need to further specify our target and projectile — since, in 
the interesting regime at high energy, the cross-section for gluon production is given 
by a rather universal formula (see Sect. 5) — it is still useful to have a more specific 
physical picture at hand, at least at intermediate stages. To that aim, we shall assume 
that the incoming hadronic systems have been produced via the high-energy evolution 
of an initial quark-antiquark {qq) dipole of size tq (within the large-iVc approximation). 
Then, the wavefunction of the projectile at rapidity yi can be described as an onium [46], 
i.e., a collection of qq dipoles which has been generated in the course of the evolution via 
1 — > 2 dipole splitting. On the other hand, the target wavefunction at rapidity 7/2 develops 
saturation effects, so its description transcends the dipole picture, even at large N,,. 

Gluon production in onium-hadron scattering can then be described as follows (see also 
Fig. 1) : a gluon radiated by any of the dipoles within the onium undergoes multiple 
scattering off the strong color fields in the hadron before eventually emerging with rapidity 
rj and transverse momentum = k. The differential cross-section for gluon production 
(or 'gluon spectrum') can be given a factorized structure [9,24-29,54]: 



which can be heuristically understood with reference to Fig. 1. The dimensionless quantity 
^('"o, f^i, yi) refers to the onium wavefunction: it represents the density of dipoles of size ri 
produced from the original dipole ro after a BFKL evolution up to rapidity yi = Y/2 — rj, 
with the initial condition n(ro,ri,0) = ro5{ri — Tq). Furthermore, agCi? ln(ri/r) can be 
roughly interpreted as the probability that an elementary dipole of size ri radiate a 
gluon-gluon (gg) dipole of size r, the dipole being effectively constructed with a gluon 
line in the direct amplitude times another one in the complex conjugate amplitude (r 
is the transverse distance between these two lines). {Tgg{r))y^ is the (average) scattering 
amplitude for the collision between the gg dipole and the hadronic target evolved up to 
rapidity 1/2 = F/2 + 77. Finally, the Bessel function Jo{kr) together with the integral over 
r implements the Fourier transform from the dipole size r to the transverse momentum k 
of the produced gluon. The presence of the derivatives w.r.t. r in Eq. (2.2), for which we 
have not given any interpretation, makes it clear that our above 'physical interpretation' 
is rather crude. A better understanding of Eq. (2.2) would require a careful inspection of 
its derivation in the literature, but this goes beyond our purposes here. 

Eq. (2.2) is often referred to as expressing "fcj_-factorization" , and it is instructive to un- 
derstand why. To that aim, it is convenient to first solve the BFKL equation for n{ro, ri, yi) 
and then insert the result into Eq. (2.2). One thus finds 




(2.2) 



drid'^k(Pb k"^ Jo 
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drJo{kr) n(ro,r, yi) 
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(2.3) 



(we have used 2Cf ~ Nc at large Nc together with ag = agNc/n), where 
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n{ro,r,y) = J^—i^-j — exp (a,|/x(7)J , (2.4) 

is essentially the Mellin representation of the BFKL solution. (^(7) is the standard BFKL 
characteristic function [48] and the contour C in the complex 7 plane is chosen as C = 
{7 = 7i + ^72 ; — cxD < 72 < 00}, with a fixed 71 in the range < 71 < 1.) Eq. (2.3) is 
finally rewritten in a form on which /c^-factorization becomes manifest : 

dcr as r dV / x^/, \ 

where 

<p{p,y) = I -^e'^-^h{ro,r,y) (2.6) 

represents (up to a normalization factor of order ctg) the unintegrated gluon distribution 
in the onium, whereas the quantity 

$(fc,y) = /dVe*'-V^(T,,(r))^ (2.7) 

can be similarly interpreted as a generalized 'unintegrated gluon distribution' for the tar- 
get, which is however sensitive to the non-linear effects in the target (related to gluon 
saturation) via the multiple scattering encoded in the dipole amplitude. It is only for suffi- 
ciently large k — much larger than the target saturation momentum — that the quantity 
{l/as)^{k,y) reduces (up to a numerical factor) to the standard, leading-twist, uninte- 
grated gluon distribution. Note also that the two distributions defined in Eqs. (2.6) and 
(2.7) have different dimensions: this is so since, whereas the onium distribution ip{p,y) 
has been already integrated over the transverse area of the onium (so this is 'unintegrated' 
only with respect to p), the target distribution ^{k, y) on the other hand is rather under- 
stood at a fixed impact parameter b within the area covered by the target (hence, this is 
more properly a gluon occupation factor in the transverse phase-space) . 

In previous applications of Eqs. (2. 2)- (2. 7), one has always treated the wavef unction of 
the dense hadronic target in the mean field approximation, that is, one has included the 
effects of saturation and multiple scattering — e.g., by computing the dipole scattering 
amplitude {Tgg{r))y with the Glauber-Mueller formula, or by solving the BK equation 
— but one has ignored the gluon-number fluctuations (although the applicability of the 
factorization (2.2) in the presence of Pomeron loops has been noticed^ by Kovchegov 
[54]). Our purpose in what follows is to go beyond such previous analyses, by including 
the effects of the Pomeron loops in the calculation of {Tgg{r))y. As we shall see, these 
effects lead to a dramatically new physical picture at sufficiently high energy. 

^ Note, however, that the analysis in Ref. [54] includes only the 'large' Pomeron loops which 
appear when the collision between the gg dipole and the target is viewed as onium-onium 
scattering in the center-of-mass (COM) frame — here, the COM of the subsystem made by the 
gg dipole and the target. Thus, this analysis is limited to not too high energies [35,46], unlike 
our present approach where the Pomeron loops are included directly in the target wavefunction. 
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(2.5) 



The Pomeron loop equations are written for dipoles made with quark-ant iquark pairs 
(rather than gluons), but this is not a problem since at high energy the scattering ampli- 
tudes for the two types of dipoles {gg and qq) are simply related to each other. Namely, 
if Sgg and Sqg denote the respective ^'-matrices, with T — 1 — S, then in the eikonal 
approximation we have: 

S,,{r) = j;^S'Jr) « S'Jr) , (2.8) 

where the last, approximate, equality holds at large Nc, and we have also assumed that 
the ^'-matrix is real, as appropriate at high energy; this in turn implies: 

Tgg{r) ^ 1 - Sgg{r) ^ 2Tgg{r) > T^^iv) . (2.9) 

Notice that, in the previous formulae, we have never indicated the dependence of the 
dipole amplitude upon the dipole impact parameter, but only upon its size. This is for 
consistency with our subsequent approximations, in which we shall systematically neglect 
the 6-dependence, i.e., we shall compute amplitudes or cross-sections at fixed impact 
parameter and assume the high-energy evolution to be quasi-local in b. 

In the next section we shaU concentrate ourselves on the qq amplitude Tgg(r) — that we 
shall denote simply as T(r) from now on — and recall some recent results concerning 
the calculation of this amplitude from the Pomeron loop equations. The consequences of 
these results for the generalized gluon distribution (2.7) and for the cross-section (2.2) 
for gluon production will be then examined in the subsequent sections. 



3 The dipole amplitude: Black spots & Diffusive scaling 

If there was not for (gluon-number) fluctuations, i.e., in the mean field approximation to 
the non-linear evolution in QCD at high energy, the dipole amplitude T(r, Y) at large Nc 
would be given by the solution to the BK equation [31, 34] — a non-linear generalization of 
the BFKL equation which is consistent with unitarity. Even in the presence of fiuctuations, 
the BK amplitude remains a reasonable approximation for the event-by-event amplitude, 
and also for the average amplitude (T(r))y at intermediate values of Y. It is therefore 
useful to begin our discussion with a brief reminder of the BK solution. 

I) Although the latter is not known exactly, its basic properties are well understood 
— thanks to a multitude of analytic and numerical studies [8,21,22,49,55-59] — and 
a piecewise analytic approximation to it can be written down. To that purpose, it is 
preferable to measure the dipole size in logarithmic units, by writing: 

T(r, Y) = r(p, Y) with P = In ^ , (3. 1) 

where ro is the unitarization scale in the target at low energy (e.g., if the target starts 
as a bare dipole at Y = 0, then is the size of that dipole). Note that large values of 
p correspond to small dipole sizes, or to large transverse momenta (p ~ Ink^) after a 
Fourier transform. 
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When viewed as a function of p for a fixed (and sufficiently large) F, the BK solution 
T(p, Y) appears as a /ront which interpolates between the unitarity (or 'black disk') limit, 
T = 1, at relatively small values of p and 'color transparency', T ~ e~^, at very large 
values of p, with the transition between the two regimes ('the front region') governed by the 
BFKL dynamics with saturation boundary conditions (see Eq. (3.2) below). The position 
of the front, conventionally defined as the value PsiX) at which T = 1/2, represents 
the unitarization scale at rapidity F, i.e., the saturation momentum QsiY) [16,60] in 
logarithmic units. Studies of the BK equation [49] or, simpler, of the BFKL equation 
supplemented with a saturation boundary condition [21,22], which reads : 

T(p, Y) = 1/2 for p = p,(y) ^ In {rlQl{Y)) , (3.2) 

reveal that, for large Y, the front propagates at constant speed: Ps{y) = ^oc^s^, with 
Ao ~ 4.88. Equivalently, the saturation momentum Qs{Y) rises exponentially with Y [60], 
with a 'saturation exponent' AqCKs. 

The most interesting region for us here is the front region at (roughly) Ps{Y) < p < 2ps{Y), 
where the scattering is weak (T <^ 1), yet it is significantly influenced by saturation, via 
the boundary condition (3.2). Within this region, the BK amplitude shows approximate 
geometric scaling, i.e., it depends upon the two kinematical variables p and Y only via 
the difference z = p — Ps{Y) = ln[l/r^Q^(F)]. More precisely, one finds [21,22] 

T(^,F) - ^e-^°^exp|-^^^^| for 1 « ^ = p - p,(F) « p,(F). (3.3) 

This expression involves the 'anomalous dimension' 70 ^ 0.63 and the BFKL 'diffusion 
coefficient' Pq = x"(7o) ~ 48.52, which are both hallmarks of the BFKL dynamics in 
the presence of saturation. Note that the last factor in Eq. (3.3), describing BFKL dif- 
fusion, depends separately upon z and Y, hence it violates geometric scaling. Still, when 
increasing Y, there is an increasing domain in z — the geometric scaling window at 
1^-2^ Zg{Y), with Zg(Y) ~ (2/5oQ;sy)^/^ — within which BFKL diffusion can be ne- 
glected, and the BK amplitude shows good geometric scaling: T ~ ze'~'^°^. On the other 
hand, for much larger value of z ^ PsO^)j one enters the regime of color transparency, 
where the BK solution reduces to the standard 'double-logarithmic approximation'. This 
is the cross-over domain towards the DGLAP dynamics and collinear factorization. 

One should stress at this point that all these remarkable predictions of the BK equation 
— the exponential growth of the saturation momentum with Y, the existence of a window 
for geometric scaling at moderately high 3> Ql{Y), and the violation of this scaling 
at even higher values of Q"^ — appear to be necessary in order to explain the HERA 
data at small-x, and the RHIC data for particle production in d+Au collisions at forward 
rapidities (cf. the discussion in the Introduction). Remarkably, these properties of the 
mean field approximation are qualitatively preserved after including the NLO corrections 
to the BFKL equation, but they suffer quantitative modifications [23], which have the 
right trend to achieve a better description of the data [15]. 

II) The basic properties of the BK solution are also preserved by fiuctuations, but only 
at the level of the event-by-event amplitude, by which we mean the amplitude obtained 
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after a single realization of the stochastic evolution. Namely, this amplitude can be given 
the following piecewise approximation [37] 



Tiz) 



1 for z < 



for l<z<L (3.4) 
e~^ for z^ L. 



where z = p — Ps(^) denotes the geometric scaling variable, as before, and we have 
kept only the dominant, exponential, behaviour in each region. As compared to the BK 
amplitude discussed previously, there are two important differences though: 

a) The front velocity A, which implicitly enters Eq. (3.4) via Ps = XagY, is considerably 
smaller than the corresponding prediction of the BK equation, Aq ~ 4.88. The difference 
Ao — A is analytically under control only in the (physically unrealistic) weak coupling limit 
[36, 37, 53]. On then obtains 

C 

^ - >'o - , 2/-, / 2^ ^ , (3.5) 

In (1/q;2) 

where the coefficient C turns out to be quite huge: C — 7r^7ox"(7o) ~ 150. Note that, when 
— ^ 0, A converges towards Aq (from the below), but only very slowly (logarithmically). 
The tendency of A to decrease when increasing is also confirmed by a study of the 
'strong noise limit', which shows that in that limit A vanishes as a power of l/a^ [61]. 

b) The front region, within which the amplitude (3.4) exhibits geometric scaling with the 
BFKL 'anomalous dimension' 70 ~ 0.63, is compact, i.e., it has a finite width L which is 
independent of Y. Specifically, L is the distance p — Ps over which the amphtude falls off 
from its saturation value T = 1 to a value of order a^, where the fluctuations wash out 
completely any mean fleld behaviour. This condition immediately yields: 

L ~ - In ^ + 0(1) . (3.6) 
70 oij 

More precisely, in the early stages of the evolution, the width of the geometric scaling 
window increases via BFKL diffusion (cf. Eq. (3.3)), until it gets stuck at its maximal 
value, equal to L. The condition -^^(l^orm) = L determines the front formation rapidity as 



agYf 



form 



2/3o7o 



(3.7) 



III) But the most spectacular effects of fluctuations appear in the calculation of expec- 
tation values, so like the average amplitude {T{p))y '■ for sufficiently large values of Y, 
this averaged quantity loses any memory of the shape of the event-by-event front (in 
particular, of geometric scaling), since it is dominated by front dispersion [38]. 

Specifically, starting with an unique initial condition at F = 0, a stochastic evolution 
generates an ensemble of fronts at rapidity Y, where the individual fronts have more or 
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Saturation / S^X 
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Diffuse^ boundary 



//,- 



Low density 



InA^cD p=lnQ^ 

Fig. 2. The diffusive saturation boundary generated by ttie evolution witfi fluctuations. 



less the same shape (cf. Eq. (3.4)) but differ from each other by a translation. That is, 
the ensemble can be described by treating the position ps of a front as a random variable, 
which on the basis of the correspondence with statistical physics can be argued to have a 
Gaussian probability distribution [37, 53, 62] : 



exp 



(ps - {ps)y 



with 



{ps) = XasY, and = 2[{pl) - (p,)^) = Dfra^Y, 



(3.8) 



(3.9) 



where the average front velocity A and the front diffusion coefficient Dfj. are analytically 
known only in the limit — >■ 0: A is then given by Eq. (3.5), whereas Dfr scales as 



V 



In^(lK) 



when a? 



0, 



(3.10) 



with a coefficient V that has been recently computed in Ref. [62]. Since the parameters 
A and Dfj. not known for realistic values of a^, in what follows we shall treat them as free 
parameters. In fact, rather than p and Y, it will be more convenient to work with the 
scaled variables z = p — (ps) and a. 

We see that, as a consequence of fluctuations, the saturation boundary become diffuse, 
within a radius cr oc y/Y around the average saturation line (ps) (see Fig. 2). This mech- 
anism should not be confused with the BFKL diffusion, which applies to the individual 
fronts, and which becomes anyway inoperative for sufficiently large rapidities Y > l^orm, 
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as we have seen. The dispersion in the saturation boundary rather means that, in the 
vicinity of the average saturation hne (ps), the hadron wavefunction is the site of wild 
fluctuations in the gluon density. When probing the hadron (say, with an external dipole) 
on a resolution scale p within the diffusive radius (|p — {ps)\ ^ c), one can either find a 
'black spot' — a gluon configuration with saturation scale Ps > p, which will completely 
absorb the projectile — , or a 'white spot' — a configuration having ps -C p, for which the 
dipole will emerge unscattered — , or any intermediate nuance of 'grey'. 

Although this picture has been obtained by working at a fixed impact parameter, it 
can be easily extended to a picture of the target in impact parameter space, under the 
assumption that the high-energy evolution is quasi-local in b (as argued in Refs. [37,38]). 
Then, the evolutions at different points b proceed (quasi)independently from each other 
and randomly generate any of the fronts which compose the statistical ensemble at Y. 
Accordingly, the average over the ensemble of fronts is tantamount to an average over 
all the impact parameters within the target area. We see that, even if one starts with 
homogeneous initial conditions at F = 0, the target becomes highly inhomogeneous at 
large F, as a consequence of (gluon-number) fluctuations in the course of the evolution. 
These inhomogeneities can be probed by a dipole projectile: A dipole with size r explores 
an area ~ around its impact parameter b. Because of the strong dispersion in the 
ensemble of fronts at large Y, a small dipole with p ^ (p^) can 'see' either 'black spots', 
or 'white', or 'grey', depending upon its impact parameter. This is illustrated in Figs. 3. a 
and b. 

The crucial feature about this situation is the fact that, for sufficiently large values of Y 
(such that cr 3> 1), the 'black spots' completely dominate the average dipole amplitudes 
[38] up to very large values of p ^ (ps), that is, including in the 'weak scattering regime' 
where (T(p))y <^ 1. When this happens, the contribution of the 'grey spots' is negligible, 
so the hadron appears to the external dipole as either black, or white (see Fig. 3.c) : mostly 
black when p ^ (Ps), but mostly white in the opposite situation where p ^ (Ps)- In fact, 
for a ^ 1 and at all the points deeply inside the diffusion radius (|p — {ps)\ <^ cr), the 
proportions of black and white are almost the same, so the average scattering amplitude 
is simply one-half ! But the dominance of the black spots actually extends up to much 
larger values of p, well outside the diffusive radius, namely so long as p — (ps) <^ cr^- 

To understand how this radically new physical picture arises, consider the average dipole 
amplitude, which in the present context is obtained as an average over the statistical 
ensemble of fronts, with weight function (3.8): 



with T{p — Ps) the single-event front in Eq. (3.4). (Higher-point correlations can be 
similarly computed.) The values of ps which dominate this convolution depend upon the 
competition between a (the width of the Gaussian distribution of the fronts) and I/7, with 
7 ~ C(l), which characterizes the exponential decay of the individual fronts. Specifically: 

i) When cr -C I/7 (the situation in the early stages of the evolution), the Gaussian 
ensemble is strongly peaked around the average front, hence the average amplitude retains 
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(a) (b) (c) 

Fig. 3. Dipole-hadron scattering in the fluctuation-dominated reginae at ^ 1. (a) a view along 
the collision axis; (b) a transverse view of the hadron, as 'seen' by a small dipole impinging at 
different impact parameters; (c) the simplified, black&white, picture of the hadron which is 
relevant for the average dipole amplitude. 

the single event profile: 

{Tip))y ^ T(p-(p,)), (3.12) 
and thus shows geometric scaling. 

ii) When a ^ I/7 (the situation at sufficiently high energy), then for all the values of 
p such that z = p — {ps) <C 70"^, the average amplitude is dominated by 'black spots', 
i.e., by configurations with pg > p, which look black (T = 1) on the resolution scale 
of the projectile. Note that, for large p > {ps), such configurations are relatively rare, 
yet they dominate the convolution (3.11) since the contributions T{p — ps) of the typical 
configurations (for which ps ~ (ps)) sxe exponentially suppressed. One then finds 

00 

{Tip))y ~ J dps Pvips) = - Erfc (^^) for - 00 < ^ < (3.13) 

p 

where the neglected terms are suppressed by, at least, one power of 1/cr and/or z/a^. 
(Note that we often use the fact that 7 ~ 0(1) to simplify the parametric estimates.) 
Eq. (3.13) involves the complementary error function. 



2 °° 



2 f 

Erfc(x) = die 



2-"^^ forx«-l 



TT \X\ 



X 



exp(— x^) 

^/7[X 



for |x| < 1 (3.14) 
for X ^ 1 . 



As anticipated, the average amplitude at high energy is approximately constant, {T{p)) 
1/2, within the large interval \z\ -C a around the average saturation line. Moreover, within 
the whole validity region of the approximation (3.13), i.e., for z -C a^, the amplitude shows 
diffusive scaling : it depends upon p (or r) and Y only via the single variable 
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Fig. 4. A 'phase-diagram' for the high-energy evolution with Pomeron loops. Shown are the 
average saturation line and the approximate boundaries of the scaling regions at large values of 
/3 ~ InQ^ (with the simplified notations (ps) = XY and = DY). Note the gradual transition 
with increasing Y from geometric scaling at intermediate energies (DY <C 1 j to diffusive scaling 
at very high energies (DY ^ 1). This transition can be more quantitatively studied on Eq. (4.3). 

^ p-{p.) ^ ln(l/r^(gg)) ^ 
a a a 

The physical regions for high-energy evolution in the kinematical plane p — Y are illus- 
trated in Fig. 4. The boundary of the geometric scaling window at large p follows from 
the discussion after Eq. (3.4). As for the diffusive scaling window, this extends up to 
Pmax — (Ps) + cr^ = (A + Z)fr)asF, cf. Eq. (3.13). A more complete expression for the aver- 
age amplitude which interpolates between these various regions will be presented in the 
next section (see Eq. (4.3)). Here, we would hke to focus on the high-energy amplitude 
(3.13), which exhibits some other remarkable properties: 

For relatively large p ^ {ps) such that 2; ^ a (with z <ti though), this function 
decreases as a Gaussian in z, rather than as an exponential. Thus, within a rather wide 
window at high Q^, whose width is increasing with Y (see Fig. 4), the scattering is weak, 
yet there is no 'twist-expansion' — i.e., no expansion in powers of ~ — for 

the scattering amplitude. 'Color transparency', meaning that (T(p)) oc e"'', is eventually 
recovered, but only for very large z ^ a"^ (cf. Eq. (4.3)). Furthermore, the amplitude 
(3.13) shows no 'Pomeron' -like growth — i.e., no power-like increase with Y at fixed p. 

Such properties might look curious, but in fact they are easy to understand in the present 
context: Both color transparency and the power-like increase with Y are properties of 
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the scattering between the projectile dipole and a dilute gluon configuration in the target. 
But as aheady explained, Eq. (3.13) rather describes the scattering off dense gluon con- 
figurations, for which the individual amplitudes have already reached the unitarity limit 
T — 1. When varying p or Y, we modify the composition of the statistical ensemble which 
is probed by the scattering, but the average amplitude changes only slowly. 

It is important to emphasize that the functional form (3.13) of the high-energy amplitude 
is extremely robust, as it follows from very few and general physical assumptions: the facts 
that the amplitude T(r) saturates the unitarity limit T = 1 for sufficiently large dipole 
sizes r and vanishes like a power of r when r — 0, and that the gluon-number fluctuations 
in the target lead to a dispersion in the ensemble of fronts, which increases with Y. By 
the same arguments, this result is universal : it is insensitive to the initial conditions at 
low energy and also to the details of the high-energy evolution. The latter matter only for 
the calculation of the parameters A and Dfr, or, more generally, of the Y-dependencies of 
the average saturation scale {ps{Y)) and of the front dispersion cr{Y). 

By using Eq. (3.13) within the dipole factorization scheme for DIS, we have been able to 
estimate (in a previous work which also included Y. Hatta and D. N. Triantafyllopoulos 
[45]) the high-energy limit of the inclusive and diffractive DIS cross-sections, and thus 
demonstrate that all the remarkable properties discussed above in relation with the dipole 
amplitude do also transmit to these cross-sections. Returning to the problem of interest 
for us here — namely, gluon production in onium-hadron collisions at forward rapidity — 
it seems natural to anticipate that these properties will show up in this context too, since 
the corresponding cross-section (2.2) is again related to the dipole amplitude. But a closer 
inspection of Eq. (2.2) may shed serious doubts on this conclusion: the 'generalized gluon 
occupation factor' defined in Eq. (2.7), and which is a building block of the cross-section 
(2.2), involves the (second order) derivative of the dipole amplitude. When taking this 
derivative on the convolution (3.11), we kill the contribution of the fronts which are at 
saturation in the statistical ensemble (the 'black spots'), because for them T = 1. Rather, 
the whole non-trivial contribution to this derivative arises from the tails of the individual 
fronts. Thus, clearly, gluon production is dominated by a different physics as compared to 
DIS, and one may wonder whether this physics is universal at high energy, and whether it 
can be computed within the present formalism. The answer to both questions turns out 
to be 'yes', as we shall argue in the next section. 

4 The generalized gluon distribution 

We start with Eq. (2.7) which is rewritten as (after performing the angular integration) 
^k,Y) = /f d--^o(^r)l|-(r|: (r,,(r))^) 

= 2 / dpJo(kroe-^/') ^(2 {T{p))y - {t\p))^), (4.1) 

where p = ln(rQ/r^) and we have used Eq. (2.9) to express the second line in term of qq 
dipole amplitudes alone. As explained at the end of the previous section, the derivative 
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dp {T{p))y is controlled by the tails of the individual fronts in the statistical ensemble. 
Hence, in order to compute this quantity, one needs to evaluate Eq. (3.11) with a complete 
profile for the single-event front, which includes a model for the tail. However, we shall 
shortly argue that the precise shape of the tail is unimportant in the kinematic range of 
interest and, moreover, the required derivatives can be simply obtained by differentiat- 
ing Eq. (3.13) ! This seems paradoxical, since Eq. (3.13) has been a priori obtained by 
neglecting the tails altogether ! This paradox is resolved by observing that the dominant 
contributions to these derivatives are actually given by the end points of the tails towards 
the saturation region (that is, by p ~ ps) 



In order to show this and also compute dp {T{p))y to the desired accuracy, we shall use 
the following, simplified, profile for a single front (compare to Eq. (3.4)) 



T{p) = 



1 for p < ps 



(4.2) 



with generic 7 ~ 0{1). A simple calculation yields (with z = p — (ps)) 
{T{p))y = I Erfc (^) + 1 exp - 7^) [2 - Erfc - 



2 



(4.3) 



The first term is the same as in Eq. (3.13) and comes from the saturation piece of the 
single front (4.2), while the second term comes from the exponential tail in Eq. (4.2). 
Incidentally, the previous results in Eqs. (3.12) and (3.13) can be easily recovered by ex- 
panding Eq. (4.3) in the corresponding limits. The corresponding expression for (^^(p))^ 
is obtained by replacing 7 — > 27 into Eq. (4.3). 

At this level, it is straightforward to compute the required derivatives by acting with 
dp = dz on Eq. (4.3). However, it is more instructive for what follows to compute at 
least the first derivative before performing the Gaussian average over p^, that is, by first 
differentiating the single-event profile (4.2) and then averaging. By using 

dpT{p) = -7e(p-p,)e-^(''-^^), (4.4) 

which implies 

p 

dp{T{p))y = -7 / dp,e-^(''-''^)Py(p,), (4.5) 
—00 

it becomes manifest that the derivative is fully controlled by the tail, as anticipated. But 
Eq. (4.5) also shows that the contributions of the fronts with p^ ^ p are exponentially 
suppressed; so, unless p is too different from (p^) (in such a way that the probability 
Py{Ps ~ p) be suppressed even stronger !), the integral in Eq. (4.5) is dominated by Ps 
close to the upper limit, within a distance I/7 from it. This yields the following estimate 



,2 

2 



dp {T{p))y ^ -Py{Ps = p) = e-^ for \z\ < a' , (4.6) 
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which can be confirmed via a direct calculation based on Eq. (4.3). Remarkably, this 
result coincides with the first derivative of the 'black spots' approximation for (T(/9))y, 
Eq. (3.13). 

Thus, as anticipated, the p-derivative of the average dipole amplitude is dominated by 
those fronts in the ensemble for which the position ps of the front is of the order of 
the dipole resolution scale. Because of that, the result (4.6) is robust and universal (i.e., 
insensitive to the details of the evolution), except possibly for its overall normalization, 
which is not fully under control in the present approximations (since this depends upon 
the precise approach of an individual front towards saturation). The same conclusions 
hold, of course, for the second-order derivative 9^, which can be computed either from 
Eq. (4.3), or directly from Eq. (3.13), with the following result: 

{T{p))y ~ e-^ for 1^1 < (7^ . (4.7) 



Since the estimates (4.6) and (4.7) are independent of 7, it is clear that they identically 
hold for the corresponding derivatives of {T'^{p))y as well. This is also in agreement with 
the fact that, within the validity range of Eq. (3.13), all the A^-body dipole amplitudes 
coincide with each other: 

(T(p)) ~ (T2(p)) ~ (T^(p)) ^ ^Erfc(^^) for a > 1 and ^ < (4.8) 

The universal functions which yield the dipole amplitude (T(p))y and its first two deriva- 
tives (appropriately rescaled by factors of a in order to exhibit diffusive scaling) in this 
high-energy regime are shown in Fig. 5. 

To complete our calculation of the generalized gluon distribution, we need to evaluate the 
Fourier transform in Eq. (4.1). We shall first give an analytic estimate, based on Eq. (4.7), 
for the dominant behaviour at high energy, and then verify this estimate, and also study a 
wider kinematical range, via a numerical calculation using the more complete expression 
(4.3) for the average dipole amplitude. 

Let us start by introducing some new notations, which will be useful in what follows: 

Z.m^. ,..„-lj = .n-5i5^-ln^ = .-Z, (4.9) 

We shall now deduce an analytic estimate for the 'occupation factor' $(A;, Y) = ^{Z, Y) 
valid in high energy regime at cr ^ 1 and for transverse momenta k within a wide interval 
around the average saturation momentum {Qs{Y)), specified by \Z\ <^ a^. By inspection 
of the integrand in Eq. (4.1), we expect the integral to be dominated by dipole sizes 
r ~ l/Zc (or z ~ Z), for which Eq. (4.7) applies: Indeed, the Fourier transform restricts 
the integration to kr ^ 1 (or z ^ Z), but the extremely small values r -C 1/A; (or very 
large z ^ Z) are suppressed by the Gaussian decay of the dipole amplitude, cf. Eq. (4.7). 

After also using the property 2 (T) — (T^) = (T), cf. Eq. (4.8), and the new notations 
introduced in Eq. (4.9), Eq. (4.1) can be rewritten as 
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Fig. 5. The average dipole amplitude in the 'black spots' approximation, Eq. (3.13), together 
with its first and second derivatives (scaled by appropriate powers of a in such a way to exhibit 
diffusive scaling) are represented as functions of the respective scaling variable t, Eq. (3.15). 



^Z,Y) = 2/dr/ Jo{e-^/') —{T{r^ + Z))^ (4.10) 

where the average ampUtude is viewed as a function of z — r) + Z — the natural variable 
in the r.h.s. of Eq. (4.7). The Bessel function Jo{x) is oscillating for large values of its 
argument x, with the amplitude of the oscillations rapidly decreasing when increasing x. 
Hence, the dominant contribution to the above integral is given by the interval between 
X — and the first zero xq of the Bessel function (i.e., from i] = r]o, with e~''°/^ = Xq, to 
77 — > oo), within which Jo{x) ~ const. This yields the following estimate 

1 f (v+z) 

$(Z, Y) ^ — j dr] (ri + Z) e ^ 

vo 

= - e > ?» - e 7^ for (7 > 1 and 1< |Z| < cr , (4.11) 

(J (7 

up to some numerical fudge factor. (The precise position of the center of the Gaussian is 
not really under control here, hence the restriction to \Z\ ^ 1.) 

Eq. (4.11) is essentially the same as the Gaussian probability distribution of the fronts, 
Eq. (3.8). This result can be physically understood as follows: ^{Z, Y) with Z = ln(A;^/ (Ql)) 
is the average 'gluon occupation factor', as obtained after averaging the corresponding 
event-by-event quantity with the Gaussian probability (3.8). In turn, the event-by-event 
'occupation factor' $event(^,^) is computed as the appropriate Fourier transform, cf. 
Eq. (2.7), of the dipole amplitude T in a single event ^ , Eq. (3.4). The result of this cal- 

^ Strictly speaking, one cannot use a piecewise interpolation, so like Eq. (3.4), in the calculation 
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iog(fcV(Q^)) 

Fig. 6. The average 'gluon occupation factor' obtained by numerical integration in Eq. (4.1 ) with 
the average dipole amphtude taken from Eq. (4.3) (with 7 = Ij for a > 2. For comparison, we 
also show the corresponding mean field curve obtained with the McLerran-Venugopalan model 
for the dipole amplitude. Note the flattening of the curve with increasing a: this demonstrates 
the rapid breakdown of geometric scaling through fluctuations when increasing the energy. 

culation is well known from previous studies based on the mean field approximation (see, 
e.g., Ref. [6]) — namely, one finds that the single-event 'occupation factor' vanishes as k"^ 
at small k <^Qs and it approaches the bremsstrahlung spectrum oc at large k ^ Qg. 
In logarithmic units, with p = ln(A;^ro), $event(p, ^) ~ e~l''~''=l is strongly peaked near 
p — Ps, and thus it behaves almost like a 5-function when averaged over ps with a Gaus- 
sian probability with a large width a ^ 1. This explains the result in Eq. (4.11). Note 
that the latter ceases to apply when \Z\ ^ cr^, i.e., for momenta k which are very much 
smaller, or very much larger, than the average saturation momentum; but in that case, 
we expect to recover the standard exponential (in p) behaviour, namely $ ~ q-\p-{ps)\_ 

This discussion allows us to put the results in Sect. 3 in a somewhat wider physical 
perspective: The results of the high-energy evolution with Pomeron loops can be visualized 
directly in terms of the (unintegrated) gluon distribution in the hadron, independently of 
the problem of the scattering. This perspective is consistent with the original derivation 
of the Pomeron loop equations [38, 39], in which the high-energy evolution was described 
as the color-glass evolution of the target wavcfunction. 

Note the contrast between the event-by-cvcnt 'occupation factor' $cvcnt(^)^)) which is 
strongly peaked at a transverse momentum k of the order of the saturation scale Qg in 



of a Fourier transform, since the discontinuities in the higher derivatives would introduce spurious 
oscillations. Rather, to that purpose one needs a smooth approximation to T(r) so like the 
McLerran-Venugopalan model [16] or the BK solution. Within such approximations, one finds 
indeed an 'occupation factor' $event(^) ^) with the behaviour discussed in the text [6, 8, 9, 30]. 
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Fig. 7. The results corresponding to a > 2 in Fig. 6 are redisplayed as a function of the diffusive 
scaling variable t = Zja. One can thus appreciate the rapid onset of the scaling behaviour when 
increasing a (or Y). The asymptotic curve is the Gaussian in the r.h.s. of Eq. (4.11). 

that event, and its expectation value ^{k,Y), which at large Y has only a mild peak at 
k ~ {Qs{Y)). This behaviour is illustrated in Fig. 6 which shows the quantity $(Z, F) 
obtained by numerical integration in Eq. (4.1) with the average dipole amplitude taken 
from Eq. (4.3). As manifest on this figure, the dispersion in the gluon distribution increases 
rapidly with a (hence, with Y), leading to the breakdown of geometric scaling. The large-cr 
behaviour is in qiialitative agreement with the analytic estimate (4.11). To more precisely 
test the latter, notice that, according to Eq. (4.11), the function (t^{Z,Y) should exhibit 
diffusive scaling at sufficiently high energy. In order to check this, we have redisplayed the 
numerical results from Fig. 6 as a function of the scaling variable t = Z/a : as manifest 
on Fig. 7, the diffusive scaling emerges indeed (and quite fast !) when increasing a. 

5 The cross-section for gluon production 

We are finally in a position to complete our original objective in this paper, namely the 
calculation of the differential cross-section (2.3) for forward gluon production in onium- 
hadron scattering at high energy. By using the expression (4.3) for the average dipole- 
target amphtude, together with the BFKL solution (2.4) for the dipole number density in 
the projectile, it is in principle possible to (numerically) evaluate the integral in Eq. (2.3) 
for arbitrary values of its external variables. However, one can always gain more insight 
via an analytic estimate, and this can be indeed obtained in the kinematical regime of 
interest for us here, namely the regime in which the Pomeron loop effects are most visible 
on this cross-section. 

Specifically, this regime is characterized by large values for the total rapidity Y and 
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also for the pseudo-rapidity r] of the produced gluon (but such that the difference yi = 
Y/2 — 1] remains relatively small, cf. Eq. (2.1)), and also by relatively large values for 
the transverse momentum k = k± of the produced gluon, in the vicinity of the average 
saturation momentum in the target {Ql{y2))- The most interesting conditions are when 
(T^ = D^6isy2 ^ 1 and \ \n{k'^ / {Ql))\ <S cr^, since in these conditions we probe the 
fluctuation-dominated regime in the target wavefunction, cf. Fig. 4. (Here and from now 
on, the quantities {Ql) and cr will always refer to the target, hence they are evaluated at 
a rapidity y2 = Y /2 + rj.) 

Under these assumptions, the convolution in Eq. (2.3) simplifies considerably, as is best 
seen on the A;_L-factorized version of this cross-section, Eq. (2.5): Given the large asym- 
metry between the evolutions of the projectile and respectively the target (recall that 
y2 ^ y\)-i the interesting values for the external momentum A;^ are much larger than the 
momentum scale Qp(yi) ~ (l/rQ)e'^"''^i which is generated in the projectile wavefunction 
by its own high-energy evolution. (Qp(yi) is the saturation scale of the projectile, and 
effectively acts as an infrared cutoff® in the integral over p in Eq. (2.5).) When this hap- 
pens, the integral in Eq. (2.5) is dominated by momenta p logarithmically distributed 
within the range Q^iyi) <^ p^ <^ k^; in the leading logarithmic approximation w.r.t. the 
transverse logarithm ln(A;^/(5p), the result can be estimated as 



da 



Oi r d'^T) 1 

drjdWb " ^^^^'^^^ / (2^^^^'^^) " ^^G,ix^,k), (5.1) 

where XiGp{xi, k"^) is the integrated gluon distribution in the projectile evaluated for a 
longitudinal momentum fraction xi — e~^^ and a transverse resolution (or momentum 
cutoff) k"^. More precisely, it can be checked that the estimate (5.1) holds under the 
following conditions: 

-cr < Z = In — — < (T^ and In — > cr . (5.2) 

The first inequality {Z ^ —a) is necessary to ensure that the contribution of the rel- 
atively large momenta p ^ k to the integral in Eq. (2.5) is indeed subleading in the 
regime where $(p, 1/2) is given by Eq. (4.11). The second inequality {Z <C cr^) ensures the 
vahdity of the approximation (4.11) for $(A;, ^2)- Finally, the third inequality is not really 
independent — it automatically follows from the previous conditions together with the 
fact that ln((Q^)/Q^) 

oc T] grows much faster with rj than a ~ — but it is nevertheless 
shown here to emphasize that the leading-twist approximation for the projectile requires 
the stronger condition \n{k'^/Qp) :$> a rather than just the usual one ln(A;^/Qp) ^ 1. 

Eq. (5.1) is recognized as the coUinear factorization of the cross-section for gluon produc- 
tion; this is vahd, as expected, at sufficiently large k. Note that this formula still involves 
the unintegrated gluon distribution of the target, ^{k, 2/2), which is moreover computed in 
a regime where the saturation effects are important. Thus, although it looks familiar, the 



® Indeed, this is the scale at which there is a turnover in (p{p,yi) from the bremsstrahlung 
spectrum oc l/p^ at large momenta to the saturation spectrum oc p"^ at low momenta. 
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collinear factorization in Eq. (5.1) emerges here in a less famihar context: the one in which 
the leading-twist approximation applies to the projectile, but not also to the target. 

Strictly speaking, the projectile distribution xiGp{xi, k"^) which enters Eq. (5.1) should be 
computed in the double-logarithmic approximation (DLA) — the approximation which 
resums both the high-energy logarithms ln(l/a;i) and the transverse ones ln(/c^ro) — , 
since it emerges from a /c_L-factorization at high energy. But this is not important here 
since, irrespective of the actual approximation which is used to compute this function 
(DGLAP or DLA), the result is anyway slowly varying with as compared to the other 
functions which appear in the r.h.s. of Eq. (5.1). 

According to Eq. (5.1), the spectrum of the produced gluon is proportional to the gluon 
occupation factor in the target $(A';,?/2) that we have already computed (in the high- 
energy regime of interest) in the previous section. To better appreciate the modifications 
introduced by Pomeron loop effects, it is preferable to consider the distribution in In k'^ : 

'^'^ «$(^,|/2), (5.3) 



dr] dln/c^ d^6 



which is simply proportional to ^{k, 1/2) ■ (We ignore the weak dependence of the projectile 
gluon distribution upon k^.) One can then distinguish various physical regimes: 

i) When k^ is extremely high, so high that the saturation effects become completely 
negligible, the gluon occupation factor in the target takes its standard perturbative form, 
i.e., the bremsstrahlung spectrum (up to the effects of the DGLAP evolution, that we 
neglect here) : 

'^'^ oc when e ^ {QD . (5.4) 



dri dln/c^ d% k"^ 

We use here the notation k"^ ^ {Ql) to emphasize that, in the fluctuation-dominated 
regime at high energy (a^ 2> 1), this leading-twist result is recovered only for momenta so 
large as compared to {Ql) that the condition Z ^ cr^ is satisfied, with Z = la(k'^ / (Ql)) . 

For lower momenta, but still much larger than the (average) saturation momentum, we 
encounter saturation effects, which can be different depending upon the energy range at 
hand (see also Fig. 4) : 

ii) For intermediate energies, such that cr^ <^ 1, the effects of fluctuations are neghgible 
and the BFKL evolution with saturation boundary conditions leads to geometric scaling 
at intermediate momenta, with an 'anomalous dimension' 70 ~ 0.63 : 

'^'^ oc f when « k' < Ql (5.5) 



dr] dlnk^ d% 



Here, we have written (Q^) = (no fluctuations !) and used Qg to denote the upper 
boundary of the geometric scaling window. We have = Q^e^^ where the exponent Zg 
grows like Zg ~ in the early stages of the evolution (cf. Eq. (3.3)), and then get stuck 
at a value 2;^ ~ L when 1/2 ^ Vform (cf. Eqs. (3.6)-(3.7)). 
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iii) For very high energies, such that o"^ 3> 1, one enters the fluctuation-dominated 
regime, where the target occupation factor $(A;,t/2) (cf. Eq. (4.11)), and hence the pro- 
duced gluon spectrum (5.3), are dominated by 'black spots' with Qg ~ k. Then, the 
average spectrum is a Gaussian in Z = hi{k'^ / {Ql)) and shows diffusive scaling (after 
multiplication by a factor of a) : 



Note the dramatic change in behaviour from the bremsstrahlung spectrum 1/A:^ oc e~^, 
which is naively expected at high k'^, to the Gaussian spectrum (5.6), which at sufficiently 
high energy should prevail up to very large values of k'^, well above {Q^ (and also well 
below it, within the range indicated in the equation above). 

6 Conclusion and perspectives 

The main conclusion of this work is that forward particle production in proton-proton 
collisions at high energy represents an ideal laboratory to search for the new physics 
expected in QCD at very high energies. Indeed, the spectrum of the gluons produced at 
forward rapidities reflects quite faithfully the unintcgrated gluon distribution in the highly 
evolved hadronic target, and thus is directly sensitive to new phenomena like gluon satu- 
ration and gluon-number fluctuations, which are expected to become important at high 
energy. The essential consequence of these phenomena is the fact that, up to relatively 
large transverse momenta — well above the average saturation momentum in the target 
wavefunction — , the cross-section for gluon production is dominated by black spots, i.e., 
rare gluon configurations with unusually large density, which are at saturation on the 
resolution scale of the produced gluon. And the distinctive signature of this behaviour, 
which should be its hallmark in the experimental results, is a new scaling law — the diffu- 
sive scaling — , which should eventually replace at sufficiently high energy the 'geometric 
scaling' behaviour currently observed in the small-x data at HERA and RHIC. 

Both geometric scaling and diffusive scaling are exceptional phenomena, in that they rep- 
resent consequences of saturation which manifest themselves at relatively large transverse 
momenta, well outside the saturation region — that is, in a region of the phase-space 
where one would naively expect the applicability of the linear, or 'leading-twist', evo- 
lution equations of perturbative QCD (DGLAP or BFKL). But although they are both 
hallmarks of saturation, these two types of scaling correspond to very different physical 
regimes, that they should unambiguously identify in the data : 

• Geometric scaling occurs at intermediate energies, i.e., energies which are large enough 
for the saturation momentum Qg to be a hard scale, yet low enough for the dispersion 
in the values of Qs (as generated via gluon-number fluctuations at low density) to 
remain negligible (o" ^ 1, cf. Fig. 4). In that regime, 'geometric scaling' is simply the 
statement that the physics remains invariant along any line parallel to the saturation 
line Ps{y) = Ay. But the dynamics when increasing p = \Rk± above Ps (i.e., when 
moving away from saturation) is still described by the leading-twist, BFKL, formalism. 




when 




(5.6) 
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• Diffusive scaling, on the other hand, represents the uhimate behaviour at sufficiently 
high energy, where the dispersion in the gluon configurations is very large (a ^ 1) and 
the high-density fluctuations dominate the expectation values within a wide window 
around the a?;era5'e saturation line {ps{Y)) — namely, so long as \p—{ps) \ <S cr^- Within 
that whole window (which includes very large momenta k± ^ (Qs) !), the leading-twist 
approximation breaks down and the cross-section for gluon production is controlled by 
the physics of saturation (similarly to the DIS cross-sections [45] ) . 

Properties of the gluon distribution produced via a high-energy evolution, the geometric 
and diffusive scaling naturally transmit (within their respective ranges of existence) to the 
scattering amplitudes for external dipoles — the most direct probes of the gluon distribu- 
tion in the target — and hence to all the processes, so hke lepton-hadron DIS or forward 
gluon production in hadron-hadron collisions, which can be given a dipole factorization. 
On the other hand, the consequences of these phenomena on the particle production at 
central rapidities — and, more generally, on the symmetric collision between two highly- 
evolved hadrons — are presently unclear, since there is no general factorization scheme 
for the scattering between two hadronic systems which are both affected by saturation. 

This problem becomes particularly acute in the high-energy regime where the Pomeron 
loop effects are important, because in that case we expect saturated gluon configurations 
to significantly contribute to scattering up to very large transverse momenta (at least, 
this is what we have seen to happen in DIS and for the forward gluon production), so the 
standard, collinear- or kj_-, factorization schemes fail to apply even at very large trans- 
verse momenta, well above the average saturation scale. It thus remains as an important 
open problem how to compute, e.g., gluon production at central rapidity in proton-proton 
scattering with Pomeron loops. It is tempting to conjecture that, within a wide kinemat- 
ical range, this process should be dominated by black-spots-on-black-spots scattering, 
but this remains to be proven. 

But further theoretical progress is also needed in the simpler case of the asymmetric, 
dilute-dense, collisions, as discussed in this paper and in Ref. [45]. The results presented 
here have been obtained via rather crude approximations, which rely on the correspon- 
dence between high-energy evolution in QCD and the reaction-diffusion process in statis- 
tical physics, but do not allow for a more detailed characterization of the QCD amplitudes. 
Although the qualitative conclusions which emerged in this way (like the dominance of 
black spots up to large transverse momenta or the diffusive scaling) are quite robust, as 
they follow from fundamental considerations, some other, equally important, properties 
— chiefly among them, the energy dependencies of the average saturation momentum and 
of its dispersion — , remain poorly understood and deserve further studies. 

To make progress, one needs more complete (possibly numerical) solutions to the Pomeron 
loop equations of Ref. [38]; this would allow one not only to determine the parameters 
A and but also to control the impact-parameter dependence of the amplitudes (to 
leading order accuracy in a.^ and 1/Nc). Furthermore, the previous experience with the 
next-to-leading-order (NLO) corrections in the context of saturation [23] shows that 
such corrections are quantitatively important for all but the asymptotically high energies. 
A NLO formahsm including Pomeron loops is therefore highly desirable, but this looks 
quite difficult to achieve. (Already the generalization of the present, LO, equations to 
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generic values of Nc meets with serious, technical and conceptual, problems, in spite of 
some recent progress [41,43,44].) There should be nevertheless possible to estimate at 
least some NLO effects, like those of the running of the coupling or of imposing energy 
conservation in the evolution. Last but not least, in view of more realistic calculations 
of the (forward) hadronic yield, it would be interesting to include some more processes 
at partonic level, so like quark pair production, and also to complete the (generalized) 
coUinear factorization, Eq. (5.1), with fragmentation functions for the produced partons. 
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